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Abstract:  A  method  and  system  of  real-time 
feedback  and  analyzing  Electroencephalogram  (EEG) 
using  short-time  Fourier  transform  (STFT)  and 
coherence  analysis  is  presented  in  this  paper.  The 
system  composition  of  hardware  and  software  are 
described.  By  the  utilization  of  this  method  and  system, 
the  EEG's  real-time  time-frequency  distribution  (TFD) 
is  obtained,  and  the  EEG  temporal  spectral  3-dimention 
array  is  displayed  on  computer's  screen  for  EEG  real¬ 
time  feedback.  The  time-frequency  analysis  (TFA)  such 
as  Wigner-Ville  distribution  (WVT)  and  wavelet 
transform  (WT),  has  been  successfully  applied  in  some 
biomedical  signals  to  detect  both  temporal  and  spectral 
features  of  biomedical  signals.  But  real-time  realization 
of  TFA  is  difficult  since  the  computation  complication 
of  WVT  and  WT.  By  using  the  STFT  algorithm,  we 
develop  a  system  for  the  EEG  feedback  and  TFA  in 
real-time.  This  system  may  be  used  for  EEG 
biofeedback,  and  possesses  the  value  of  clinical 
practice. 
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1.  INTRODUCTION 

Biofeedback  is  a  training  technique  in  which 
people  are  taught  to  alter  brain  activity,  blood  pressure, 
muscle  tension,  heart  rate  and  other  bodily  functions 
that  are  not  normally  controlled  voluntarily,  and  to 
improve  their  health  and  performance  by  using  signals 
from  their  own  bodies. 

EEG  biofeedback  involves  helping  a  person  learn 
how  to  modify  his  or  her  brainwave  activity  to  improve 
attention,  reduce  impulsivity,  and  to  control  hyperactive 
behaviors.  It  is  a  painless,  non-invasive  treatment 
approach  that  allows  the  individual  to  gain  information 


about  his  or  her  brainwave  activity  through  the  use  of 
computerized  biofeedback  equipment  and  use  that 
infonnation  to  produce  changes  in  brainwave  activity. 
EEG  Biofeedback  is  an  effective  treatment  for 
headaches,  attention  deficit  hyperactive  disorder 
(ADD/ADHD),  epilepsy,  sleep  onset  disorders  and 
insomnia.  More  than  700  groups  worldwide  are  using 
EEG  biofeedback  (neurofeedback)  for  treatment  of 
ADD/ADHD. 

Biomedical  signals,  such  as  electrocardiogram 
(ECG),  electroencephalogram  (EEG),  event-related 
potentials  (ERP),  electromyogram  (EMG),  etc.,  are 
fundamental  observations  for  analyzing  body 
functions  and  for  diagnosing  diseases,  and  are  useful 
for  both  physiological  research  and  medical 
applications.  Information  provided  by  those  signals  are 
generally  time-varying,  nonstationary,  sometimes 
transient,  and  usually  corrupted  by  noise. 

Time-frequency  analysis  (TFA)  has  been 
successfully  applied  in  some  biomedical  signals.  The 
main  transforms  of  Cohen's  class  allow  signal 
representation  simultaneously  in  time  and  frequency 
domain.  There  are  various  methods  to  transform  signals 
into  Time-frequency  distributions  (TFD).  Recent 
research  in  Time-frequency  analysis  have  concentrated 
on  Short-time  Fourier  transform  (STFT),  Wigner-Ville 
distribution  (WVD),  Exponential  distribution  (ED), 
Discrete  Gabor  spectrogram  (DGS),  Wavelet  transform 
(WT),  and  Wavelet  networks  (WN).  Cross-term 
interference  occurred  when  using  the  WVD  algorithm. 
This  interference  seriously  disrupted  the  measurement 
of  such  biosignals,  despite  the  fact  that  the  WVD 
algorithm  showed  better  time  resolution.  Choi  and 
Williams  introduced  a  distribution  called  the 
exponential  distribution,  in  which  cross  terms  can  be 
controlled  by  a  single  parameter. 

Wavelet  transform  has  better  resolution  comparing 
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to  short-time  Fourier  transform  and  has  not  the  cross¬ 
term  interference  comparing  to  the  Wigner-Ville 
distribution.  It  has  better  features  in  time  -frequency 
domain.  Wavelet  networks,  a  combination  of  wavelet 
transform  and  artificial  neural  networks,  represent 
signals  as  linear  wavelet  nodes  which  characterized  by 
time-frequency  features  related  to  the  wavelet 
transform. 

But  real-time  realization  of  TFA  is  difficult  since 
the  computation  complication  of  WVT  and  WT.  By 
using  the  STFT  algorithm,  we  develop  a  system  for  the 
EEG  TFA  and  feedback  in  real-time.  This  system  may 
be  used  for  EEG  biofeedback,  and  possesses  the  value 
of  clinical  practice. 

2.  MATHEMATIC  AT.  METHODS 

Short-time  Fourier  transform  (STFT)  is  one  of  the 
Time-frequency  analyses.  The  STFT  transform  of  one 
channel  EEG  signals  h(n),(n=0,l,...,N-l),  is  given  by 

H(k)  =  J^w(n)h(n)e-2mklN 

n= 0 

Where  w(n)  is  Kaiser  Window. 

The  alpha  index  (the  EEG  alpha  rhythm  relative 
density)  is  defined  by 

I  II  =P  II  /P  EEG 

I  a  =P  a  /P EEG 

I  fj  =P  p  /Peeg 

here  P  a  and  PEEG  are  the  EEG  alpha  rhythm  power  and 
total  power  respectively. 

For  the  computation  of  coherence  spectra,  consider 
two  channels  of  EEG  data,  x(n)  and  y(n),  (n=0,l,...N-l). 
By  STFT  analysis,  one  obtains: 

X(k)=Ax(k)+jBx(k) 

Y(k)=Ay(k)+jBy(k) 

From  these  amplitudes  A(k)  and  B(k),  the  cross¬ 
spectra  C  and  Q  are  determined  from  the  relations 
Cxy(k)=Ax(k)Ay(k)+Bx(k)By(k) 
Qxy(k)=Bx(k)Ay(k)-By(k)Ax(k) 

To  improve  the  statistical  properties  of  these  cross- 
spectral  estimates,  it  is  in  general  necessary  to  smooth 
them  in  some  appropriate  fashion.  Denoting  this 
smoothing  process  by  boldface  type,  the  coherence 


function  at  frequency  kf/N,  is  defined  by 


clXV+QjXk) 

cjvcjk) 


We  may  also  define  a  total  coherence  function  for 
the  band  [kl,k2]  by 


£[ clm+Qim 
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In  order  to  draw  COSPAR  on  a  printer,  the 
coherence  function  should  be  dealt  with  as  following 


rjk)= 


J  Yxy  (k)  -  h, 

jo, 


7jk)>h 

7jk)<h 


Here  h  is  the  threshold  volume. 

We  define  the  EEG  dominant  frequency  fmax  and 
alpha  index  I„  as  two  EEG  feedback  object  functions 
and  they  are 

I  c  =I  c  (F3)+I  a  (F4)+  T  F3  F4(kl  ,k2) 

Where  I  a  (F3)  and  I„(F4)  are  alpha  index  of  EEG's 
lead  from  F3  and  F4  areas,  is  the  alpha  band  coherence 
between  F3  and  F4. 


fmax=fmax(F3)+fmax(F4) 

Where  fmax(F3)  is  the  dominant  frequency  of  F3 
and  fmax(F4)  the  dominant  frequency  of  F4.  In  the 
feedback  strategy,  two  criteria  are  used.  The  first  is  to 
increase  the  index  Ia,  and  the  second  is  to  decrease  the 
index  fmax. 

The  automatic  artifact  rejection  is  the  key  aspect  for 
the  EEG  closed-loop  feedback  control  and  it  determines 
the  robust  and  stability  of  the  feedback  system.  By 
comparing  the  fmax  with  the  baseline  volume,  the  system 
decides  whether  the  artifacts  are  serious  and  automatic 
rejection  should  be  carried  out. 


3.  SYSTEM  CONFIGURATION 


The  EEG  real-time  feedback  system  consists  of 
EEG  amplifier  ,  EEG  pre-processor,  A-D  converter, 
IBM  PC/AT  microcomputer,  earphone.  Printer,  CRT. 
Printer  and  CRT  are  used  to  output  EEG  feedback 
control  results.  EEG  preprocessor  is  a  notch  filter 
inserted  at  50Hz. 
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Fig2.  The  EEG  COSPAR  (h=0.8) 

System  software  consists  of  3  parts  that  are 
initiation,  EEG  signal  acquisition  and  analysis,  result 
output.  Firstly,  and  the  Kaiser  window  w(n)  is 
calculated.  Then,  the  EEG  signals  are  digitized  on-line 
to  12  bits  at  64  samples  per  channel  in  a  A-D  converter 
(PC- 1232)  connected  to  a  IBM  PC/AT  microcomputer. 
Simultaneously  with  the  data  acquisition,  eight 
channels  of  EEG  are  processed  by  adding  Kaiser 
window,  and  short-time  Fourier  transform.  By  use  of 
EEG  power  spectrum,  EEG  alpha  rhythm  relative 
power  (alpha  index),  EEG  main  frequency,  and  total 
coherence  in  alpha  band  are  calculated,  the  EEG 
temporal  spectral  3 -dimention  array  is  displayed  on 
computer's  screen  for  EEG  real-time  feedback  (Figl). 
On  the  other  hand,  according  to  those  EEG  indexes,  the 
computer  will  determined  the  strategy  of  feedback,  and 
control  the  earphone  to  send  out  the  auditory  signals 
for  EEG  feedback. 

Finally,  the  system  outputs  the  EEG  index  dynamic 
curve  and  the  EEG  coherence  spectral  array  (COSPAR) 
(Fig2)  on  the  printer  and  CRT. 

4.  RESULTS 

Subjects  were  trained  with  the  EEG  feedback 


system.  They  were  seated  in  a  comfortable  chair  within 
a  sound-controlled  room.  With  eye  slight  closing, 
subjects  concentrated  on  the  warning  sound  from  the 
earphone,  and  tried  to  decrease  the  appearance  of  the 
sound  by  subjective  effort.  Records,  with  linked  ears 
for  reference  electrode,  were  taken  from  electrodes 
F3/F4  of  the  10-20  system.  The  ground  was  to  Fz. 

In  the  early  days  of  training,  subjects  were  not  able 
to  effectively  regulate  themselves  EEG's  activity  due  to 
lack  of  experience.  After  several  sessions  of  training, 
they  gradually  mastered  the  ability  to  control  EEG's 
activity. 

The  experimental  results  of  EEG  feedback  training 
suggested  that  EEG's  activity  was  able  to  be  effectively 
regulated  after  several  sessions  of  training,  and  the 
effect  of  TM  practice  was  received  by  the  utilization  of 
the  EEG  feedback  system. 
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